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A physical scenario accounting for the existence of size-limited submicrometric domains in cell 
membranes is proposed. It is based on the numerical investigation of the counterpart, in lipidic 
membranes where proteins are diffusing, of the recently discovered cluster phases in colloidal sus- 
pensions. I demonstrate that the interactions between proteins, namely short-range attraction and 
longer-range repulsion, make possible the existence of stable small clusters. The consequences are 
explored in terms of membrane organization and diffusion properties. The connection with lipid 
rafts is discussed and the apparent protein diffusion coefficient as a function of their concentration 
is analyzed. 
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I. INTRODUCTION 

The question of membrane functional organization is a 
key issue in modern cell biology [H, 0, S IIS S, 0, H, ©]• 
The central problematics is to establish the relationship 
between (dynamical) organization and functions of the 
different constituents of the membranes. In this context, 
many microscopy techniques are implemented to have ac- 
cess, with the highest possible spatial and temporal res- 
olutions, to the distribution and dynamics of membrane 
proteins and lipids, in particular to their diffusive prop- 
erties in connection with their crowded environnement. 

Using these techniques, there exist a large variety of 
situations, in live cells [E II 0, i, 0, [H (3 or in model 
membranes [U, H!)] , where membrane proteins are found 
in oligomers or small clusters. This co-localization is sup- 
posed to facilitate encounters of partners in different pro- 
cesseSjSuch as complexes involved in signal transduction 
0, 0, ll- Lipid rafts @, HI HH are usually invoked to 
account for co-localization. These membrane submicro- 
metric domains enriched in certain lipids {e.g. choles- 
terol and sphingolipids) are supposed to recruit proteins 
having a higher affinity for their composition. They 
are believed to ensue from a lipidic micro-phase sepa- 
ration 0, HJ, HH . However, a consensus has not yet been 
reached to explain why the separation process stops and 
domains remain size-limited [151 ] . And understanding the 
relationship between these structural patterns and diffu- 
sion of their constituents remains topical P, 0, [f| 0, [ll| . 

At the same time, there is an increasing interest in 
colloid science for systems presenting a cluster phase. It 
is the fruit of a competition between a short-range at- 
traction {e.g. a depletion interaction) which favors clus- 
terization, and a longer-range repulsion {e.g. electro- 
static) which prevents a complete phase separation, be- 
cause when clusters grow, their repulsion also grows and 
the repulsion barrier cannot be passed by thermal acti- 
vation anymore. The result is an equilibrium phase with 
small clusters of concentration-dependent size. It has 
been suggested [l6| that the existence of cluster phases 
is not a singular behavior in a specific system. They in- 
deed occur in different physical systems {e.g. colloids, 
star polymers, proteins) and for a large variety of inter- 
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have already been experimentally observed or simulated 
in 2D 0,[2i,[2l[25i. 

I propose, as an alternative paradigm, to re-interpret 
the aggregation of membrane proteins in terms of two- 
dimensional cluster phases, because membrane proteins, 
like colloids, interact with energies of the order of mag- 
nitude of the physio logical thermal energy ksT {T « 
310 K) [1, IMS III III- I discuss in this framework the 
mechanism driving the formation of so-called rafts. I pro- 
pose that proteins spontaneously congregating in small 
clusters instead of constituting a "gas" of independently 
diffusing inclusions, they promote the formation of nano- 
domains in membranes. The sub-micrometric, limited 
size of these nano-domains now appears naturally in this 
scenario, as a result of the competition between attrac- 
tion and repulsion. By contrast to the lipid raft scenario 
where domains result from a lipidic micro-phase separa- 
tion and then recruit specific proteins, the present mecha- 
nism proposes that domain formation is mainly driven by 
protein interactions. Recent experiments in live cells sup- 
port this point of view that protein-protein interactions 
are necessary to induce clustering, whereas lipids alone 
are not sufficient in the system studied [IH . I also relate 
the protein diffusive properties as a function of their con- 
centration to the limited size of clusters, itself depending 
on concentration. I demonstrate that the mean long- 
term diffusion coefficient D of proteins decreases when 
the mean cluster size grows: D oc l/{n) where (n) is 
the mean cluster size (its number of proteins). I also an- 
ticipate that (n) grows with the protein concentration 4> 
[l3|]. Thus D also depends on tf> [ll|. From this point 
of view, contact is made with prior experiments, where 
D was found to decrease like 1/0 [3(| or that were inter- 
preted by appealing to such a behavior [j|. I propose a 
simple scenario leading to this law D ~ l/(f>, requiring 
that (n) cx (f>. Previous attempts to address this issue, 
based on low <fi expansions of D, are not adapted to catch 
the physical mechanisms capable of accounting for this 
behavior (3lj . 

To give substance to and to support this cluster-phase 
scenario, I first discuss available experimental data. In 
Ref. [HI], proteoliposomes of egg lecithin and bacteri- 
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orhodopsin (BR) are observed by freeze-fracture electron 
microscopy. BR is found to aggregate in small clusters 
(called "particles"), the mean size of which, (n), depends 
on BR concentration cf>. In Table I (column 5), it is found 
that (n) ex <j> for l/<j> S [90,300] (in mol./mol. lipid 
to protein ratio), in complete agreement with my argu- 
ments. In addition, the Fig. 4 of this reference shows un- 
ambiguously that this behavior comes from the bimodal 
character of the cluster size distribution: the balance be- 
tween monomers (n = 1) and multimers (n > 1) make 
possible the cluster density to be nearly constant while 
the density of protein grows. I will confirm below that it 
is exactly what is observed in simulations. Moreover, the 
diffusion constant D of BR is precisely found to decay 
like l/<p in the same interval in Ref. [301 ] . as anticipated 
above. At the end of this article, I shall propose experi- 
ments to validate definitely this scenario. 



II. SHORT AND LONGER-RANGE 
INTERACTIONS 

There exist several short-range attractive forces be- 
tween proteins embedded in membranes, each with a 
range of a few nanometers and a binding energy of or- 
der ksT. They first of all consist of a depletion in- 
teraction due to the 2D osmotic pressure of lipids on 
proteins, which tends to bring them closer when they 
are about a nanometer apart [261 ]. There also exist hy- 
drophobic mismatch interactions between proteins, the 
hydrophobic core of which does not match the width of 
the membrane [2|. The energy cost of the subsequent 
membrane deformation increases with the distance be- 
tween two identical proteins, thus resulting in an attrac- 
tive force. The energy scale is of order fceT [2fij |. In 
membranes with several lipid species, another scenario 
leads to attractive forces: proteins recruit in their neigh- 
borhood lipids which match best their hydrophobic core. 
The closer the proteins, the more energetically favorable 
the configuration. Binding energies are also of order ksT 
or larger [2, 29] . A protein-driven mechanism for domain 
formation invoking such forces has been proposed 0, H^] , 
but the limited domain size due to additional repulsive 
forces has not been discussed in this context. 

Membrane inclusions are also affected by longer-range 
repulsive forces. Electrostatic repulsion between like 
charged proteins is usually considered as negligible be- 
cause it is screened beyond a few nanometers in physio- 
logical conditions: at physiological ionic strength I v ~ 
0.1 M, the Debye screening length is of the order of 
1 nm [321 ]. Only proteins with (unreasonable) charges 
of several hundreds of elementary charges can give a re- 
pulsion of a fraction of ksT at 10 nm. By contrast, there 
exist specific repulsions due to the deformation that pro- 
teins impose on the membrane when they are not strictly 
speaking cylindrical inclusions but conical ones or periph- 
eral proteins (27j . For example, using the formulae of this 
reference for transmembrane proteins with a moderate 



contact angle of 10°, one finds that the repulsive energy 
barrier at 10 nm is 0.10 ksT for a typical bending rigid- 
ity k = 100/csT. For example, rhodopsin has a contact 
angle larger than 10° [33| . 

Thus the ingredients for the existence of cluster phases 
are present in assemblies of membrane proteins and clus- 
ter phases should generically exist in cell membranes. Be- 
low, I shall take a typical binding energy of —4 ksT (l7j 
and an energy barrier of 0.1 fe^T at about 10 nm. 



III. CLUSTER DIFFUSION 

First I study the diffusive properties of particles in clus- 
ter phases. I consider an isolated cluster of n proteins, 
modeled as an assembly of interacting Langevin particles 
with "bare" diffusion coefficient Do (the diffusion coef- 
ficient at vanishing concentration). The center of mass 
of the assembly diffuses with a coefficient Dq/u, because 
the clusters considered are not rigid entities but loosely 
bound, fluctuating ones in which the proteins diffuse [5| 
(this property will be confirmed below) . If clusters inter- 
act weakly because they are sufficiently far away (Fig.[T]), 
the long-term diffusion coefficient of each protein of the 
cluster is also equal to Do/n [5(. If clusters contain (n) 
proteins on average (counting a monomer as a cluster 
with n = 1), then the mean long-term diffusion coeffi- 
cient D = Do/(n): if the system contains N proteins 

c— 1 x ' 

where N c — N/ (n) denotes the number of clusters, 
because a cluster c contains n(c) proteins that diffuse 
each with a diffusion coefficient Do/n(c). Thus D — 
D N C /N = D /(n). If R is the average cluster radius, 
then D oc 1/R 2 . Such an experimental behavior of D 
with cluster size has already been observed [3], but has 
been left unexplained. 

Now, in 3D, it has been shown analytically (l8l. l2ll. [231] . 
and measured experimentally [H, E3, El] , that (n) is pro- 
portional to the particle concentration 4>. My purpose 
here is not to demonstrate such a relation, but simply 
to remark its validity in a wide range of situations in 
3D and to anticipate its equivalent in 2D. In addition 
to the numerical evidence presented below, further cal- 
culations, appealing for example to the theory of micel- 
lization [341] . will be necessary to confirm this last point. 
They go beyond the scope of the present numerical pa- 
per. If (n) oc 4>, then the effective diffusion coefficient 
of proteins in a cluster phase is inversely proportional to 
their concentration: 



D = Const. /(j). (2) 
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FIG. 1: Two snapshots of the cluster phase at <f> = 0.1 
(TV = 200 proteins, box side a = 0.25 /im, periodic boundary 
conditions); The time delay between both snapshots is 0.5 ms. 
Clusters diffuse slowly and appear non-rigid, deformable and 
fluctuating. 



IV. MONTE CARLO SIMULATIONS 

A priori, these mechanisms should be valid only at 
low <j), where separate clusters interact weakly and dif- 
fuse independently. To test further their relevance in 2D 
and at higher 0, I have performed Monte Carlo simu- 
lations of systems of N = 100 particles (N = 200 for 
the highest density considered tfi = 0.1), interacting via 
pairwise potentials displaying a hard-core repulsion, a 
short-range attraction and a longer-range repulsion [2J|. 
I have chosen physically and biologically relevant param- 
eters as justified above. I have observed a strong robust- 
ness of the cluster phase with respect to the potential 
shape. I have tested repulsive terms decaying linearly, 
algebraically (like 1/r 2 or 1/r 4 (2?]]) or exponentially [2l[ 
with the distance r between molecules, as well as attrac- 
tive ones varying exponentially or linearly with r. In all 
cases, a cluster phase exists at equilibrium (i.e. after 
very long runs) for wide ranges of parameters. Clus- 
ters co-exist with a gas of monomers, of density depend- 
ing on potential and (weakly) on density 4> (as observed 
experimentally in Ref. .13]). What is important is not 
the precise shape of the potential but the existence of 
a short-range attraction of a few fceT and of a longer- 
range, weaker repulsion extending on a range larger than 
the typical cluster diameter. Therefore I have focused on 
a potential shape already studied in detail [2l|, [22|, [24[ : 



U(r) = -e a exp(-7 a r) + e r exp(-7 r r). 



(3) 



The parameters are chosen so that, as required above, 
the binding energy between two proteins is —4 fc^T and 
the energy barrier is 0.1 fcsT. The following values ful- 
fill this requirement: e a — 32 fc^T and e r = 0.3 fc^X 1 ; 
l/7 a = 2 nm and l/7 r = 16 ran. In spite of the high 
value of e , the binding energy is low because the attrac- 
tive part is cut at r — do due to the hard core repul- 
sion (Inset of Fig. [5]). This hard-core diameter is cho- 
sen as d$ = 4 nm, the typical diameter of a protein of 



average molecular weight [1J|. The proteins are given 
a "bare" diffusion coefficient D = 1 fim. 2 /s [3(J l36j : 
at each Monte Carlo step (MCStep), a randomly cho- 
sen protein attempts to move a distance Sr forward in a 
randomly chosen direction; Here 5r = 1 A <C do; With 
this 5r, the acceptation rate of MCSteps is larger than 
60 %, even at the highest densities considered; The time 
step between two MCSteps is St — <5r 2 /(4£> ) = 2.5 ns. 
A Monte Carlo sweep corresponds to TV MCSteps. The 
simulation time is chosen so that error bars on D and 
(n) are smaller than 10 % (more than 10 7 sweeps, i.e. 
30 ms of real time). The protein average long-term diffu- 
sion coefficient D is measured at different concentrations 
<\> = Nd^/a 2 (a is the size of the box with periodic con- 
ditions in which the proteins diffuse) . By long-term, it is 
meant at time-scales larger than the time needed to dif- 
fuse inside the clusters (typically 0.1 ms). The measures 
are performed after a long equilibration period. To be 
sure that equilibrium has indeed been reached, I sim- 
ulate two systems with initial configuration chosen as 
(a) a random one where proteins later coalesce to form 
clusters; (b) a completely condensed state where all pro- 
teins belong to the same big cluster which later splits 
into smaller ones and gas. Equilibrium is considered to 
have been reached when both systems (a) and (b) are 
qualitatively identical (same number of multimeric clus- 
ters). The Monte Carlo time needed is generally shorter 
than 10 7 sweeps (however, see [HI). Note that the time 
needed to reach equilibrium in (a) is rather long because 
after a transient period where small clusters appear via 
a bimodal-like decomposition, larger clusters are formed 
by evaporation of the smaller ones. Evaporation is the 
result of escape of single proteins from the clusters, one 
after the other. The energy barrier to evaporate a single 
protein being of several fc^T, this is a slow process [35j . 



V. RESULTS AND DISCUSSION 

As 4> was increased, D decreased dramatically, as ex- 
pected (Fig. [2]). One observes that D « Const./ <p over 
nearly two decades, thus confirming the hypothesis that 
cluster phases can account for this behavior. This law 
had been observed 25 years ago [3(3| , without receiving a 
full explanation, apart from arguments invoking "crowd- 
ing effects" or "aggregation" reminiscent of the clustering 
mechanism discussed here, but unable to predict quan- 
titatively the dependence of D on <fi [30]. The diffusion 
coefficient at typical cell membrane protein concentra- 
tion (</; ~ 0.1) appears to be reduced by a factor larger 
than 10, which is the decrease of D observed in cell mem- 
branes as compared to the same diffusion coefficient Do 
in model membranes at low <f> 0, [33] • One can also see in 
Fig. [2] that (n) is proportional to <fi at low concentration, 
as expected. Interestingly, at higher </>, the data for (n) 
still follow this law for this choice of potential. In addi- 
tion, the relation D ~ Do / (n) holds on all the range of 
concentrations studied. This demonstrates that the in- 
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FIG. 2: The diffusion coefficient D/Do (diamonds) and the 
inverse mean cluster size 1 / (n) (circles) as a function of the 
inverse density l/<j>, in log-log coordinates, for the potential 
U discussed in the text (see [35| for <f> = 0.1). The full line has 
slope 1, for comparison. Inset: the potential U(r) (in units of 
feT; r in nanometers). 



FIG. 3: Numerical cluster-size distributions, for three differ- 
ent concentrations: 4> = 0.0031 (black), <f> — 0.0125 (red) and 
4> = 0.05 (gray). For 4> — 0.0031, multimers are not visible 
with the scale used because they are very scarce. The frac- 
tions P(n) of clusters of size n are in the same units as cj>: 
number of clusters per unit surface do- 



teractions between clusters are negligible and that they 
diffuse independently at the time scale considered. 

As it was observed experimentally (Fig. 4 of Ref. [IH), 
the cluster distributions obtained in simulations are bi- 
modal on all the range of concentrations where (n) oc <j>. 
A gas of monomers (n — 1) coexists with large mul- 
timers (n > 1), the distribution of which is Gaussian 
around a typical size n* . With the parameters chosen 
here, there are virtually no small multimers (dimcrs, 
trimers, etc...). These distributions are illustrated in 
Fig. [3] at different concentrations 0. The relation (n) oc <f> 
then comes from a subtle balance between monomers and 
multimers: as 4> increases, the density of monomers is es- 
sentially constant while large clusters capture additional 
proteins. 

Once formed, clusters appear to be non-rigid, de- 
formable, fluctuating (as illustrated on the two consec- 
utive snapshots in Fig. [TJ they can be seen as liquid 
droplets, that reorganize rapidly; proteins diffuse inside 
the clusters with measured short-term diffusion coeffi- 
cients larger than 0.01 /im 2 /s) and long-lived, as antici- 
pated above. By long-lived, it is meant that clusters are 
stable at the time-scale of the simulations (about 30 ms) . 
However, some proteins constantly leave the clusters (via 
the evaporation process discussed above) , diffuse freely in 
the gas phase and are captured later by another cluster. 
Rare events of clusters being disintegrated or nucleating 
spontaneously in the gas have been observed. Therefore 
clusters are certainly not stable at long time scales (sec- 
onds or minutes) and are only transitory. 

I have also observed in simulations that small mod- 
ulations of the parameters e or 7 in U(r), can lead to 
segregation [l|, [T^, [26[: if two (or more) groups of pro- 
teins (A's and B's, which are not necessarily identical 



in the same group) are present in the simulation and if 
A-A and B-B interactions are slightly favored as com- 
pared to A-B ones, then the proteins segregate. There 
are A-rich and B-rich clusters because even though it 
is entropically unfavorable, it is energetically favorable. 
For instance, an A-B binding energy 10 % smaller than 
equal A-A and B-B ones suffices to ensure segregation. 
At the biological level, this implies that groups of pro- 
teins which show a slight tendency to associate, because 
of their physico-chemical properties, would segregate in 
distinct clusters (see [Hj]). This mechanism might play 
an important role in sorting together proteins that must 
congregate to perform biological functions [l], This 
point will be quantified in future investigations [34j |. 

The existence of cluster phases in live cells would mean 
that proteins spontaneously congregate in small clusters 
of a few entities or a few tens of entities in the plasma 
membrane instead of constituting a "gas" of indepen- 
dently diffusing inclusions, thus promoting nano-domains 
in plasmic membranes. By contrast to the lipid raft sce- 
nario where domains result from a lipidic micro-phase 
separation and then recruit specific proteins, the present 
mechanism proposes that domain formation is mainly 
driven by protein interactions (even though lipids do play 
an important role in the effective forces). Note that this 
scenario does not exclude a concomitant recruitment by 
protein clusters of specific lipids having a higher affin- 
ity for those proteins (and which participate in the ef- 
fective attractions 29]), thus reconciling my hypothe- 
sis with the observation of detergent-resistant membrane 
fractions [lij]. This mechanism, by constraining an in- 
creasing fraction of lipids to diffuse slowly (with clus- 
ters) as <j) increases, could also explain why the diffusion 
constant of lipids decreases significantly when the con- 
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centration of proteins increases (30j . 

In spite of the evidence provided above, the existence 
of protein clusters has to be confirmed definitively at 
the experimental level. Even though such clusters have 
already been observed by freeze-fracture electron mi- 
croscopy [3, [HI , their existence must be explored by dif- 
ferent techniques in a wider range of situations, in cell 
and model membranes. Near-field scanning microscopy 
is an ideal tool because it is able to identify individual 
proteins after immobilizing the membrane onto an ade- 
quate substrate. Counting numbers of proteins in clus- 
ters is then in principle possible |sj. Such experiments 
would be able to investigate the dependency of cluster 
numbers (n) on protein concentration, as well as the cor- 
relation between (n) and D. Recently, high-frequency 
single-particle tracking has demonstrated that proteins 
are confined in nano-domains in the plasma membrane 
of live cells @ , of typical diameter a few tens of nanome- 
ters. An appealing hypothesis is that these nano-domains 
are the clusters under consideration. For proteins several 
nanometers apart, the previous size would correspond to 
clusters containing a few tens of proteins, in agreement 



with the previous simulations. The confirmation of this 
hypothesis would provide additional evidence of cluster 
phases in live cells. 

I have proposed in this paper a paradigm leading to 
the formation of size-limited nano-domains in cell and 
model membranes. This scenario, based on reasonable 
hypotheses about the energy and length scales in biolog- 
ical membranes, sheds light on several so-far open issues 
in cell biology: (i) it provides a mechanism for the lim- 
ited size of nano-domains; (ii) it gives a complete qualita- 
tive justification for former experiments in model mem- 
branes [H,[33|; (iii) it proposes a simple explanation lead- 
ing to the proportionality law D cc (j). If cluster phases 
were to be experimentally confirmed in model membranes 
and in live cells, it would mean that, by physical mecha- 
nisms, proteins generically gather in small assemblies in 
biological membranes, thus shedding new light on mem- 
brane functional processes. 
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